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We analytically investigate the time-optimal unitary evolution of entanglement between indirectly 
coupled qubits in a trilinear Ising chain with unequal interaction couplings. The intermediate qubit 
is controlled via a local magnetic field. We find the time-optimal unitary evolution law, and we 
quantify residual entanglement via the two-tangle between the indirectly coupled qubits, for all 
possible sets of initial pure quantum states of a tripartite system. The integrals of the motion of 
the brachistochrone are determined by fixing the minimal time at which the residual entanglement 
is maximized. Entanglement plays a role for W and GHZ initial quantum states, and for the 
bi-separable initial state in which the indirectly coupled qubits have a nonzero value of the 2-tangle. 
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I. INTRODUCTION 

The concept of entanglement is one of the key fea- 
tures distinguishing the quantum world from the classi- 
cal world, as it captures those correlations which cannot 
have a classical origin [1 . It is one of the fundamen- 
tal resources in quantum information and computation 
theory, teleportation, superdense coding, quantum cryp- 
tography, black hole physics etc Several measures of en- 
tanglement have been proposed for multipartite quan- 
tum systems in different contexts, for a review see, e.g., 
[2]- [5]. On the other hand, the speed of evolution for 
quantum systems is an important concept, not only for 
determining the theoretical limits at which quantum in- 
formation can travel [6] -[7], but also for the practical 
task of building quantum computers capable of perform- 
ing fast quantum algorithms before the ubiquitous and 
disruptive decoherence effects come into play. The im- 
portance of the connection between quantum entangle- 
ment, the speed of evolution of quantum systems and 
dynamical optimization problems has been discussed in 
[8]- [12]. Furthermore, quantum optimal control is also a 
fundamental subject, both theoretically and experimen- 
tally, in quantum computing and information (see, e.g., 
[13 ). In particular, time-optimal quantum computation, 
where the cost to be optimized is the time to achieve a 
given quantum evolution is relevant for the design of fast 
elementary gates and provides a more physical ground 
to describe the complexity of quantum algorithms. A 
theoretical framework for the quantum brachistochrone 
(QB) was introduced in [14 . The QB [ST is based on 
a variational principle enforcing the time-optimal evo- 
lution of a quantum system whose Hamiltonian is sub- 
ject to a set of constraints (e.g., a finite energy, certain 
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qubit interactions are forbidden) and defines a bound- 
ary value problem with fixed initial and final quantum 
states (or unitary transformations). The QB has been 
studied for quantum state evolution in the case of pure 
[14 and mixed states [15 , for the optimal realization of 
unitary transformations up to a given target quantum 
gate [16], and for the more realistic situation where the 
target can be reached within a finite, tolerable error (a 
fixed fidelity) [17 . Some authors HI]- [19] started the 
study of the role of quantum entanglement during the 
QB evolution of multipartite distinguishable systems in 
a pure quantum state, finding that the entanglement is 
pivotal to the QB evolution if at least two subsystems 
actively evolve. Efficient generation of random multipar- 
tite entangled states has been also analyzed with the aid 
of time-optimal unitary operations [20 . More recently, 
the authors of [21 found that genuine tripartite entan- 
glement is necessary during the QB evolution of a set 
of three qubits in the pure state, except for the case in 
which less than three qubits attend evolution. 

In this paper we discuss the role played by the entan- 
glement in the QB formalism by considering the time- 
optimal unitary evolution of a tripartite system with a 
trilinear Ising interaction Hamiltonian with a fixed en- 
ergy available and where the intermediate party can be 
controlled via a local, time dependent magnetic field. We 
point out that the time-optimal synthesis of interactions 
between qubits which are indirectly coupled via an in- 
termediate qubit is a typical scenario in a wide array of 
promising (scalable) experimental realizations of quan- 
tum information processing (see, e.g., [22 -[26 ). More- 
over, the QB formalism naturally allows for the situation 
in which local coherent controls are assumed to time con- 
suming, contrary to the requirements of zero time cost for 
local controls typical of the standard time-optimal quan- 
tum control methods [2T, We concentrate our attention 
on the behavior of the residual entanglement between the 
indirectly coupled qubits at the end of the trilinear chain, 
as expressed by the so called 2-tangle [28 . We also con- 
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sider all the possible initial quantum pure states for the 
tripartite system, i.e. completely separable, bi-separable, 
and with true tripartite entanglement (bipartite for W 
and tripartite for GHZ states). We then let the 2-tangle 
evolve along the general time-optimal quantum trajec- 
tory defined by the QB action principle, and we fix the 
integrals of the motion by imposing that entanglement 
reaches its maximum value in the shortest possible time. 
The minimal time for reaching such a maximum of the 
2-tangle is a function of the ratio between the interaction 
couplings of the Ising Hamiltonian. 

The paper is organized as follows. In Section II we 
review the main features of the QB formalism for the 
time-optimal synthesis of unitary quantum evolutions. In 
Section III we summarize the QB solution for the prob- 
lem of a three-linear qubit system subject to an Ising 
interaction with unequal couplings and a local control 
on the intermediate qubit, when a finite energy is avail- 
able. In Section IV we define the measure for the bi- 
partite entanglement between the two indirectly coupled 
qubits of the chain and we introduce the main formulas 
for the 2-tangle and 3-tangle for all the possible sets of 
initial quantum states. Section V is devoted to the study 
of the entanglement evolution for these initial quantum 
states, and we define the optimal times and the analyti- 
cal form of the quantum evolutions for which the 2-tangle 
between the indirectly coupled qubits most rapidly reach 
its maximum value. Finally, Section IV is devoted to the 
summary and discussion of our results. 

II. TIME-OPTIMAL UNITARY EVOLUTION 

The goal is to determine the time-optimal way to gen- 
erate the unitary evolution up to a certain Uf (mod- 
ulo physically irrelevant overall phases) by controlling an 
Hamiltonian H{t) obeying the Schrodinger equation. We 
assume that H is controllable and that only a finite en- 
ergy is available in the experiment. This time-optimality 
problem may be formulated using the action 

S{U,H;a,A,Xj):= f ds[Na ^ Ls ^ Lc] , (1) 

JO 

Ls:={K,i%U^ -aH), (2) 
Lc'.= aY,^,P{H), (3) 
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where {A^B) := Ti^A"^ B) and the Hermitian operator 
K{s) and the real functions \j{s) are Lagrange multipli- 
ers. The quantity a is the time cost, relating the param- 
eter time s and the physical time t via t \= J a{s)ds. 
Variation of Ls by A gives the Schrodinger equation: 

= HU, or U{t) = Te-'^o^dt^ 

where T is the time ordered product. Variation of Lc by 
\j leads to the constraints for H: 

fi{H) = Q. (5) 



In particular, the finite energy condition for a system of 
logA^ qubits reads: 

m):=\[T,{H^)-NLo^]=0, (6) 

where is a constant. 

From the variation of S with respect to H we get: 

A = Aoff + 5:A,^, (7) 

while from the variation of 5* by a we obtain the normal- 
ization condition: 

Tr(i^A) = N. (8) 

Finally, variation of 5 by /7, use of eq. ([7| and some el- 
ementary algebra give the quantum brachistochrone equa- 
tion: 

.f = [^,A], (9) 

The quantum brachistochrone together with the con- 
straints define a boundary-value problem for the evo- 
lution of the unitary operator U{t) with fixed initial 
{U{t = 0) = 1, where 1 is the identity matrix) and final 
conditions {U{t = T) = Uf^ where T is the optimal time 
duration necessary to achieve the target gate Uf). It can 
be solved together with the constraint functions fj{H) 
to obtain i^opt(^)- Then one integrates the Schrodinger 
equation Q with /7(0) = 1 to get Uopt{t) and finally the 
integration constants in i^opt(^) can be fixed, e.g., by im- 
posing that Uopt (T) equals a target U f modulo a global 
(physically irrelevant) phase. 

III. QB FOR A TRILINEAR ISING CHAIN 

We now apply the general formalism of the QB to the 
case of a physical system of three qubits (labeled by a 
superscript a G {1,2,3}) interacting via an Ising Hamil- 
tonian and where the intermediate qubit is subject to a 
local and controllable magnetic field Bi{t) {i = x^y^z): 

Hit) := Ji2a\al + Jasa^cr^ + B{t) ■ a\ (10) 

where we have defined cr^^cr| (g) <jj (g) 1, cr|<j| := 

1 (g) (Ji (g) CTj, cr? := 1 (g) (g) 1 and cr^ are the Pauli 
operators. Introducing the ratio between the Ising in- 
teraction couplings K := J23/J12 and rescaling time as 
r := Ji2t, the energy as Co := UJ / J12 and the magnetic 

field as Bopt{t) — Ji2^opt('?")7 it can be shown [16] -[17] 
that the QB is solved by the following time-optimal mag- 
netic field: 

f Bo cos 0{r) \ 
4pt(r)= ( 5osin^(r) 1 , (11) 



3 



where 0{r) := flr^Oo and 1^, Oq, Bq and are integra- 
tion constants. The magnetic field is precessing around 
the 2;-axis with the frequency Ct. Furthermore, the energy 

constraint exphcitly reads = cl;^ — (1 + K'^) -= ujI^ so 
that we can reparameterize: 



Bq := cos (j); B^ := u^k sin ( 



(12) 



where and (j) G [0,27r]. 

One can then solve the Schrodinger equation Q and 
find the time-optimal evolution operator [16]-[T7]: 

C/opt(r) = e-^^l,3(r)|0)2(0| 

- ie''^-''--^BoS},^{T)\l)^{0\ + II.C., (13) 

where we have introduced the diagonal operators acting 
in the Hilbert space of qubits 1 and 3 

A]^{t) := Diag[ai(T),a2(T),a3(T),a4(T)], (14) 
Sj^ir) := Diag[si(T),S2(T),S3(T),S4(T)]. (15) 

The latter operators depend upon the functions of time: 

sin(a;^r) 



Si{r) :-- 

UJi 

Ci{r) := cos(cjir), 

ai{r) := ^(r) + i biSi{r), 

and on the constants: 

UJi := cjfe^cos^ (j) + 6f , 
1 



(16) 

(17) 
(18) 



bj := sin ( 



+ {Si2-5i3Kl-K)-^- 
where i = 1, 2, 3, 4 and 6ij is the Kronecker symbol. 



(19) 



IV. ENTANGLEMENT IN THE 1-3 
SUBSYSTEM 

For a tripartite quantum system in a pure state |?/^) = 
EI=o^*K) expanded in the basis {|0) := |000), |1) := 
|001),|2) := |010),|3) := |011), |4) := |100), |5) := 
|101),|6) := |110),|7) := |111)} the entanglement be- 
tween two of its subsystems, e.g. 1 and 3, may be quan- 
tified by the 2-tangle [5]. [28]: 



ri3 := 2[Det[pi] - Det[p2] + Det[p3] - |HypDet(a) 



(20) 



where pi (i = 1, 2, 3) are the reduced density matrices 



Pi 

P2 

P3 



= TV23(|^)W|), 



(21) 



and HypDet(a) is Galley's hyperdeterminant for the ma- 
trix of the coefficients a^s (see eq. (46) in the Appendix). 
Furthermore, one may also define the 3-t angle: 



Ti23 := |HypDet(a) 



(22) 



which describes the amount of tripartite entanglement 
between all the spins. 

We now consider the time optimal evolution of an ar- 
bitrary initial pure state |'^(0)) driven by the evolution 
operator (13), and we study the behavior of the bipar- 



tite entanglement between the indirectly coupled qubits 
1 and 3. We are interested, in particular, about deter- 
mining the optimal time and the integration constants 
Bq^Bz-, 1^, ^0 for which the 2-tangle ri3 is maximized. 

As it is well known, pure states of a tripartite quantum 
system can be classified into equivalence classes under lo- 
cal operations and classical communication (LOCC), dis- 
tinguishing the degree of entanglement between its sub- 
systems [29 . In particular, we will consider the repre- 
sentatives of each of these classes as a possible initial 
quantum state: 

a) initial completely separable states, for which all 2- 
tangles and the 3-tangle vanish: 



|V(0)) = \'^s) ■■= |000); 



(23) 



b) initial bi-separable states, with bipartite entangle- 
ment, for which one of the 2-t angles is nonzero while the 
other 2-tangles and the 3-tangle are zero: 

|^(0)> = I^Bsi) := ^(1001) + 1010)), (24) 

|^(0)> = l^BSi) ■■= ^(I001> + 1100)), (25) 

1^(0)) = l^BSi) ■■= ^(I010> + 1100)); (26) 

c) initial Werner (W) states, with full bipartite entan- 
glement, for which all 2-tangles vanish while the 3-tangle 
is nonzero: 

1^(0)) = \i>w) ■■= ^(1001) + 1010) + |100)); (27) 

d) initial Greenberger-Horne-Zeilinger (GHZ) states 
with maximal tripartite entanglement, for which all 2- 
tangles and the 3-tangle are nonzero: 

IV'(O)) = li^GHz) ■■= ^(1000) + 1111)). (28) 



The time-optimal evolution operator t/opt('?"), eq. (13), 
drives the initial states according to 



|V'(r)) = C/opt(r)|V(0)) = 5]a,(r)|i). 



(29) 



i=0 



For each of the possible classes of initial states, we com- 
pute the nonzero time-dependent amplitudes, which are 
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explicitly given by formulas (51)-(56) of the Appendix. 



Then it is easy to check, substituting these formulas into 
(|46|) and (ISOl), and then into ^ and that both the 



^isi'T) and the ri23(r) tangles are always zero during the 
whole time-optimal evolution of a 3-qubit state which is 
initially fully separable or a bi-separable state of type 
Bl and B3 (i.e., when initially there is no entanglement 
between the indirectly coupled qubits 1 and 3). 

On the other hand, when there is some initial entan- 
glement between qubits 1 and 3, then we have a non 
trivial time-optimal evolution of the ris{r) and ri23(r) 
tangles. For example, when the initial state belongs to 
the class B2 of bi-separable states, from eqs. (20), (|22|), 



but simple algebra we find that "optimal" time at which 
ri3^2 i'^) ^^^^ reaches its maximum depends on the value 
of the available energy cu'^ and on the ratio K between 
the couplings in the Ising Hamiltonian. In more details, 
we have that: 



V3 



4 1-K 



(35) 



evolutions: 



when 1 < < 29/16 and \K\ 
uj^ > 29/16 and Ki_ < K < K2- 

where we have defined Ki± 
K2± := (1/4) [13/4 ± V'3(c2;2 -29/16)]. xhe duration 



< or when 

or K2+ < K < 
and 



1 



(46), (50) and (53) and we get the following time-optimal (43) of the optimal quantum evolution for the bi-partite 



(30) 



ri23B2(^) = ^l«2S3 - a3S2p = 1 - ri3^^{r). (31) 

The case of an initial W state is similar and we find that 

riswi^) = (4/9)n3B2(^) and ri23w(r) = (4/9)ri23s2 
Finally, for the class of fragile, fully entangled GHZ 
states the time-optimal evolution of the tangles is found 
to be given by: 



(32) 



entanglement is minimal for a ratio of the Ising couplings 
K Ki±, and maximal for a ratio of the Ising couplings 
K K2±. Furthermore, in this case we obtain the op- 
timal magnitudes and frequency for the magnetic field 



\BzB2 

Qb2 = 2 



^/3 



\K-l\, 



3 3 21' 



7 8 4Q 

K-1±xIlo^- -if2 +°K- — 
3 3 21 



(36) 
(37) 
(38) 



Instead, we have that: 



V. TIME-OPTIMAL EVOLUTION OF 
ENTANGLEMENT 

The behavior of the time-optimal evolution of the tan- 
gles ri3,ri23, which witness the entanglement between 
qubits 1 and 3, can be studied in more details by ex- 
plicit use of the formulas (16), (18) and (19) into equa- 
tions ([30| 



(33). As the tangles are always related by 
1 — '?"i23(t), we limit ourselves to the study of 
the ri3 tangle. 

For the case of an initial bi-separable state of class B2 
(and similarly for initial states of the class W), from (18) 
and ([30l we obtain: 



^13b2{^) = [C2C3 +CJ^(C0SV + ^2^3)5253]^ 

+ ^fc(^252C3 - b3S3C2f. (34) 

We now proceed by further optimizing the tangle 
ri3^2 i'^) as a function of the time r and of the unknown 
constants of the motion for the time-optimal quantum 
trajectory, i.e. ^o,^^,!^, for a given ratio of the cou- 
plings K and for fixed energy cu. In particular, we look 
for the time at which the tangle first reaches its max- 
imum value ri3|max = 1 by imposing that the partial 
derivatives of the tangle with respect to r and the other 
constants of the motion (expressed in terms of (/) and Cl) 
are zero and that the determinant of the Hessian ma- 
trix at this point is negative. After a long and tedious 



2K 



(39) 



when c2;2 > 29/16 and K2- < K < K2+. Within this 
range of allowed energies and the duration of the op- 
timal quantum evolution for the bi-partite entanglement 
is minimal for a ratio of the Ising couplings K i^2-, 
and maximal for a ratio of the Ising couplings K i^2+ • 
In this case, the optimal values for the magnetic field are: 



JB2 
ZB2 



\B,^A=VL = Q. 



(40) 
(41) 



In Figs. 1-2 we plot the 2-tangle and the 3-tangle as 
a function of time when the quantum system is in the 
initial bi-separable state B2 and for the value of en- 
ergy uj = a/6- For this value of energy, we can com- 
pute Ki± = ±2.24 and i^^2+ ^ 1-70, K2- - -0.007. 
To exemplify the values of the couplings, we two models 
proposed in |25j for nuclear magnetic resonance (NMR) 
experiments, i.e.: a) the H — N — H chain in the molecule 
of ethanamide, for which J12 = J23 — 88.05 Hz and 
therefore i^^ = 1; b) the P — F — H chain in the 
molecule of diethylfiuoromethylphosphonate, for which 
J12 - 46 Hz, J23 - 73.1 Hz and therefore K = 1.59. 
For both K = 1 and K = 1.59 we are in the sit- 
uation depicted by formula ([39|, and the the law of 
quantum evolution of the tangle is explicitly given by 
ri3s.{r) = 1 - 4[(1 - Kfojl/iw^ - 2f] sin'^Wco^ - 2Kt]. 
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We notice that the case of equal Ising couphngs, i.e. case 
a) with K = 1, residual entanglement between the in- 
directly coupled qubits is constant and always maximal, 
equal to one, while the tripartite entanglement is always 
zero. For the case of unequal couplings, instead, the bi- 
partite entanglement has a periodic behavior, starting 
from the maximum equal to one at r = 0, reaching a 
minimum of approximately 0.57 at r = rising 
again to the maximum of one at r = r*^2 o^- 

Let us now turn to the case of an initial state belonging 
to the GHZ class. From (ulsl) and (1321) and we obtain: 



^ISghz (^) = COS^ 0[(C1S4 - C45i)^ 

-^ulib^-bifisis^f]. (42) 

In this case, the "optimal" time at which Ti3^^^(r) first 
reaches its maximum is found to be: 



V2 TT 



(43) 



provided that o)^ > 3/2 and K_ < K < where 
we have defined K± := (l/2)[-l ± V2uj^ - 3]. For the 
GHZ initial states, the duration of the optimal quantum 
evolution for the bi-partite entanglement is minimal for 
a ratio of the Ising couplings K and maximal for 

a ratio of the Ising couplings K K_ (and it diverges 
if cj^ = 2, and K K- = — 1). The associated optimal 
magnitudes and frequency for the magnetic field in the 
GHZ case are: 



available, we have found the time-optimal unitary evolu- 
tion law for the Ising Hamiltonian plus the local control 
and we substituted it in the formula for the 2-t angle. The 
initial boundary condition for the QB is chosen among all 
possible sets of tripartite quantum pure states, i.e. fully 
separable, bi-separable, Werner and GHZ states. The 
integrals of the motion in the QB are determined impos- 
ing that the 2-tangle reaches its maximum in the shortest 
time possible, which we call r*. Entanglement is found to 
have a non trivial role during the time-optimal unitary 
evolutions of W and GHZ initial quantum states, and 
of the bi-separable initial state in which the indirectly 
coupled qubits have a nonzero value of the 2-tangle. The 
optimal time also sets the time-scale for the duration 
of the significant role of the entanglement, and it is a 
function of the ratio K between the interaction couplings 
in the Ising Hamiltonian and of the energy available in 
the experiment. The monogamy of entanglement shows 
neatly in the anti-correlation of the tripartite entangle- 
ment, quantified by the 3-t angle, with the bipartite en- 
tanglement, quantified by the 2-tangle. It is interesting 
to extend the analysis presented here to the case when 
coherent control is possible on all the qubits in the chain, 
to longer chains of qubits, and to study the QB evolution 
of the truly non classical correlations via their measure, 
quantum discord [30]. 



\Bo,„,\ = \l + K\ 



GHZ] 



(44) 

y^uj^-2{K^^K^l). (45) 
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VII. APPENDIX 



For a tripartite system 123 in a pure state |?/^) = 
X]i=o^*K) expanded in the basis {|0) := |000), |1) := 
|001),|2) := |010),|3) := |011),|4) := |100), |5) := 
|101),|6) := |110),|7) := |111)}, Cailey's hyperdetermi- 
nant for the matrix of the coefficients a.s is defined as 



HypDet(a) := [(aoay)^ + (ai^e)^ + (<^2tt5)^ + (<^3«4)^] 
— 2[{aoar + aiaQ){a2a^ + 0.3(^4) 
-h aoaiaear + a2asa4a^] 
-h ^{aoasa^aQ + aia2a4a7). (46) 



More explicitly, since 
Det(/9i) 



|aip(|a4|2- 
\a2?{W?- 
\az?{\a,\^- 



-\a^? + \a,\^) 
lael' + KI') 
\a^? + \a7?) 
lasl' + lael^) 



+ 01070305 + 02070305), 



Det(p2) = |aoP(|o3|2 
+ |aiP(|a2|'- 
+ |a4|^(|a4|^- 

+ |a5|'(|04|2- 



03 



«7 



|a3|' + |06|2) 



and 



+ aiaja%^a^ + a4a7a5ag), 



Det(p3) = |aoP(|o3|2 + |o5p + |o7p) 

+ |a2p(|ai|' + |a5|' + |«7|') 

+ |a4p(|ai|' + |a3|' + |a7|') 

+ |a6|'(|ai|' + |a3|' + |a5|') 
— 2Re(aoa3a^a2 + aQa^a[al 

+ a2a7a%aQ + a/^ajaloQ), 



we obtain: 



(47) 



(48) 



(49) 



Det(pi) - Det(/92)+Det(/93) = 2(|oo|2|o5|2 + \ai\'^\ai\^ 

+ \a2?\ai? + \a^?W\^) 

— 2Re[aoa7a^ae + a2a^a%^aX 

— {aoaj - aiaG){a2al - a^al)] 

— 4Re(aoa5aJa4 + a2a7a'^aQ). (50) 



Given the time-optimal unitary evolution (29), we can 



compute the following nonzero time-dependent ampli- 
tudes (all modulo exp[— zl^r/2]): 

a) representative of the class of fully separable initial 
states, IV^s-): 



{cLs)o{r) = ai, 
{as)2{r) = -iBoe'^si; 



(51) 



bl) representative of the class of bi-separable initial 



states, \iIjbsi)' 



(am)o(r) = '^^^'''''i^ 
a* 

(aBi)i(r) = -|, 
(aBi)2(r) = 

(«Bi)3(r) = "^^^'^^2; 



(52) 



b2) representative of the class of bi-separable initial 
states, \1pB2)' 



(aB2)i(r) = 

\/2 



(aB2)3('r) = -i-^e'"s2, 



(aB2)4(r) = 



(aB2)6(T) = -i-^e"'s3; 



V2 



(53) 



b3) representative of the class of bi-separable initial 
states, \ipB3)' 



^0 ^-iOo 



51, 



(aB3)2(r) = -^^1. 

a* 

(aB3)4(r) = 

(«b3)6(t) = 



(54) 



c) representative of the class of Werner initial states, 
{aw)o{r) = -i^e-'^'si, 



{ciw)i{r) = 

{ciw)3{r) = -i^e'^S2, 
a* 

{aw)4{r) = 



(55) 
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d) representative of the class of GHZ initial states, 

\'^GHZ)- 



{ciGHz)o{r) 



( \ ( \ ie 



V2' 
.Bo 
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FIG. 1: The 2-t angle (red curve) and the 3-t angle (blue curve) 
as a function of time for the B2 initial states and the coupling 
ratio K — 1. 
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FIG. 3: The 2-tangle (red curve) and the 3-tangle (blue curve) 
as a function of time for the GHZ initial states and the cou- 
pling ratio K = 1. 



